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ABSTRACT 

^ ' Context. X-ray cavities and extended radio sources ('cocoons') surrounding active galactic nuclei (AGN) have been detected by 

the Chandra X-ray mission and radio interferometers. A joint analysis of X-ray and radio maps suggests that pressure values of 
non- thermal radio-emitting particles derived from the radio maps are not sufficient to inflate the X-ray cavities. We propose using 
the Sunyaev-Zel'dovich (SZ) effect, whose intensity strongly depends on the pressure, to find the hitherto undetected, dynamically- 
dominant component in the radio cocoons. 

Aims. Numerical simulations show that plasma with a high temperature (10' - 10'" K) is a good candidate for inflating the AGN 
, cocoons. To constrain the population of high energy electrons inside AGN cocoons that is predicted by numerical simulations, we 

HH ' study different methods for maximizing the contribution of such energetic electrons to the SZ effect. 

MH ' Methods. Our calculations of intensity maps of the SZ effect include relativistic corrections and utilize both analytic models and 

numerical 2-D simulations. 

^ ^, Results. We demonstrate that the spectral function at a frequency of 217 GHz has an absolute maximum at a temperature higher 

I . than 10' K, therefore the measurement of the SZ effect at this frequency is a powerful tool for potentially revealing the dynamically- 

O ' dominant component inside AGN jet-driven radio cocoons. A new method is proposed for excluding the contribution from the low 

energy, non-relativistic electrons to the SZ effect by means of observations at two frequencies. We show how one may correct for 
C/3 , a possible contribution from the kinematic SZ effect. The intensity maps of the SZ effect are calculated for the self-similar Sedov 

solution, and application of a predicted ring-like structure on the SZ map at a frequency of 217 GHz is proposed to determine the 
energy released during the active jet stage. The SZ intensity map for an AGN cocoon in a distant elliptical is calculated using a 2-D 
numerical simulation and including relativistic corrections to the SZ effect. We show the intensity spectrum of the SZ effect is flat at 
^ , high frequencies if gas temperature is as high as k\,T^ = 500 keV. 

, Key words. Radiative transfer - relativistic processes - cosmology: cosmic microwave background - galaxies: cluster: intracluster 

■ medium - galaxies: ISM 

in 

■ 1. Introduction A relativistically correct formalism based on the probability 
, distribution of the photon frequency shift after scattering was 

O . Inverse Compton scattering of the cosmic microwave back- ^^^^^ Wright (1979) to describe the Comptonization pro- 

T-H . ground radiation (CMB) by free electrons in intervening ionized ^^^^ ^f soft photons by mildly relativistic electrons. The Wright 

■ gas allows us to study the physical state of the ionized gas (for a formalism (as is the Kompaneets approximation) is based on 
review, see Sunyaev & Zel'dovich 1980). The Comptonization assumptions: (1) The Thomson cross-section is applicable; 
process of a soft photon spectrum by a Maxwellian distnbu- (2) the photons are sufficiently soft that the energy of the 

H , tion of non-relativistic electrons is described by the Kompaneets scattered photons is less than that of the electrons. The vaUd- 

>^ , equation (Kompaneets 1957). An analytical solution of the ^^^^^^ of the Wright formalism was studied by Loeb et al. 

Kompaneets equation was derived by Zel'dovich & Sunyaev ^ ^gg^^ ^^o have shown that this formahsm is vahd for many as- 

(1969) and apphed to the distortion of the cosmic microwave trophysical environments. The relativistically correct Wright for- 

background radiation by Thompson scattering in clusters of j^^jj^j^ ^^^^ calculate the SZ effect in galaxy clusters by 

galaxies (the Sunyaev-Zel'dovich effect, hereinafter the SZ ef- ^^^^ ^^^^ors (for a review, see Rephaeli 1995 and Birkinshaw 

feet). The SZ effect has proven to be an important tool for cos- 1999) 

mology and the study of clusters of galaxies (for a review, see Although the SZ effect is observed in clusters of galaxies, 
Birkinshaw 1999). It measures the pressure of an electron popu- j^ere are other environments in which this effect should be ob- 
lation integrated along the line of sight and does not depend on gervable by the next generation of telescopes with high sensitiv- 
spectral features of the underlying electron distributions, as long g^oj angular resolution (such as the ALMA telescope 
as they are non-relativistic. ^ Since a significant fraction of the energy released by active 

galactic nuclei (AGN) goes into heating of the neighbouring 
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plasma by shocks, distortion of the CMB due to the SZ effect 
by shock-heated plasma may be produced. 

According to the standard evolutionary scenario for strong 
double radio sources (Scheuer 1974; Blandford & Rees 1974), 
jets boring through the intergalactic medium (IGM) are not in di- 
rect contact with the undisturbed IGM, but rather are enveloped 
in a cocoon consisting of shocked jet material and shocked IGM. 
Scheuer (1974) showed that the cocoon around a pair of super- 
sonic, low-density (when compared to the ambient IGM) jets 
acts as a "wastebasket" for most of the energy deposited by the 
jets. Begelman & Cioffi (1989) argued that the cocoons in many 
observed sources have not yet had time to reach pressure balance 
with the ambient medium. 

The Comptonization of the CMB by the hot plasma around 
strong radio sources was considered in the framework of the non- 
relativistic Kompaneets approximation by Yamada et al. (1999), 
Platania et al. (2002) and Cavaliere & Lapi (2006). However, 
numerical simulations by Antonuccio-Delogu & Silk (2008) and 
Tortora et al. (2009) show that the temperatures can become very 
high within the cocoons {T ^ 10**- 10" K), and therefore consid- 
eration of the relativistically correct SZ formalism is necessary. 
Kino et al. (2007) found that thermal electron temperatures of 
the AGN cocoons are predicted to be of the order of an MeV and 
determined only by the bulk Lorentz factor of the jet. Therefore, 
the relativistic Wright formalism must be applied to correctly 
derive the distortion of the CMB by AGN cocoons. The tem- 
peratures in AGN cocoons are orders of magnitude higher than 
those of clusters of galaxies (T^i ~ 10^ - 10** K). 

Recently, numerous X-ray cavities in the intra-cluster 
medium have been detected by the Chandra X-ray observatory 
(e.g. McNamara et al. 2000; McNamara et al. 2005). They of- 
ten coincide with the radio lobes of the central radio galaxy. 
The non-thermal pressure derived from the equipartition con- 
dition for the energy of synchrotron-radiating non-thermal elec- 
trons and magnetic fields is a factor of ten smaller than the pres- 
sures required to inflate the bubbles (e.g. Ito et al. 2008). This 
implies that most of the energy in the cocoon is carried by an 
invisible component such as high energy thermal electrons (e.g. 
Ito et al. 2008). Observations of the SZ effect were proposed 
by Pfrommer et al. (2005) to probe the inferred dynamically- 
dominant component of plasma bubbles associated with X-ray 
cavities. Pfrommer et al. (2005) studied examples of several 
different physical scenarios concerning the composition of the 
plasma bubble which is as a whole in pressure equilibrium with 
the ambient ICM, and derived the SZ effect due to electrons in 
the plasma bubbles using the Wright formalism. However, co- 
coons are expected to be over-pressured with respect to the am- 
bient IGM (Begelman & Ciofli 1989), and therefore the SZ ef- 
fect derived under the assumption of pressure balance gives only 
a lower limit on the true SZ effect produced in the cocoons. 

To constrain a population of high energy electrons inside 
AGN cocoons predicted by numerical simulations, we study in 
this paper the induced CMB distortion as a function of gas tem- 
perature. Taking into account the derived temperature depen- 
dences, we consider different methods for maximizing the con- 
tribution of such energetic electrons. By using the self-similar 
Sedov solution and the results of our numerical simulations, we 
calculate the relativistically correct SZ intensity maps. 

The plan of the paper is as follows. We consider the SZ effect 
by mildly relativistic electrons in the framework of the Wright 
formalism in Sect. 2 to find a convenient method for observ- 
ing an electron component with high temperatures (T » 10** 
K) derived from numerical simulations of AGN cocoons. The 
Sedov stage of the AGN cocoon evolution is considered in Sect. 



3, where intensity maps of the SZ effect are calculated. The SZ 
effect produced by high energy electrons residing in a simulated 
AGN cocoon is studied in Sect. 4. We show in Sect. 5 that the 
CMB spectral distortion is flat in the broad high frequency range 
if gas temperatures are high. We calculate an intensity map of 
the SZ effect by means of the dynamical scaling of cocoon in 
Appendix A. 



2. Generalized spectral functions. 

In this section, the spectral properties of the SZ effect are studied 
within the framework of the Kompaneets and Wright formalisms 
to determine the contribution of the high temperature electron 
component. 



2.1. A generalized spectral function 

The CMB spectral distortion produced by the SZ effect in the 
Kompaneets approximation is (Zel'dovichA Sunyaev 1969): 



AI„r(x) = /og(jc)3'gas 



(1) 



where x - hv/kbTcmb, Iq - '^■{khT cmhf' I {hcf' , and the spectral 
function g{x) is given by: 



X exp(x) / exp(x) + 1 



(exp(jc) - 1)^ \ exp(x) - 1 



(2) 



The subscript '«r' denotes the fact that Eq. ([TJ was obtained in 
the non-relativistic limit. The Comptonization parameter jgas is 
given by 



^gas 



meC- J 



(3) 



where the line-of-sight integral extends from the last scattering 
surface of the CMB radiation to the observer at redshift z=0, Tg 
is the electron temperature, «gas is the number density of the gas, 
and cTj is the Thomson cross-section, the electron mass, c 
the speed of light, the Boltzmann constant and h the Planck 
constant. 

The CMB spectral distortion in the relativistically corrected 
formalism can be written in an analogous form to that given by 
(EnBlin & Kaiser 2000; Pfrommer et al. 2005), and is 



A/(x) = /,) 



(4) 



Note that in our notation, the spectral function g{x) is changed 
to the generalized spectral function G{x, Tg) which depends ex- 
plicitly on the electron temperature. 

The relativistic spectral function G{x, T^) derived in the 
framework of the Wright formalism is given by 



G{x, Te) 



Pi(s,T,)( x^exp(-3s) 



©(Te) \exp(jcexp(-i)) - 1 exp(x) - 1 



\ds{5) 



where ©(Te) = kbTg/nisC-, fi(s, Te) is the distribution of fre- 
quency shifts for single scattering (Wright 1979; Birkinshaw 
1999). 

The following limiting cases hold: 

1. if Te «; 10^ K, the values of relativistic G{x, Tg) and non- 
relativistic g{x) spectral functions are equal; 
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2. if » lO'^K, the value of the relativistic spectral function 
is written as 



G{x,T,) ^ — 



nifC 



exp(x) - 1 ^bT'e 
where i{x) - -jc-^/(exp(x) - 1) is the Planck spectrum. 



(6) 



Note that in the papers by EnBlin & Kaiser (2000) and Pfrommer 
et al. (2005), a different form for the relativistically correct 
Wright formalism is proposed. These authors have chosen to av- 
erage over the gas temperature along the line of sight. The two 
forms coincide if the astrophysical medium is homogeneous. 
However, if the astrophysical medium is inhomogeneous, use of 
the form given by Eq. (01 is more convenient since this uses local 
temperature values. 

In order to study the properties of the generalized spectral 
function G(x, Te), we calculate this function at frequencies 90 
GHz and 217 GHz at which the SZ effect will be observed 
by the ALMA telescope. These frequencies are also observ- 
able by other telescopes (like the South Pole Telescope and the 
Atacama Cosmology Telescope) and are motivated by the fol- 
lowing reasons: 90 GHz is close to the frequency of the SiO (sil- 
icon monoxide) line, and 217 GHz corresponds to the crossover 
frequency of the non-relativistic SZ effect. 

Pfrommer et al. (2005) simulated a Green Bank Telescope 
(GBT) observation of the central regions of the Perseus and 
Abell 2052 clusters where X-ray cavities have been detected 
by the Chandra observatory, and computed the frequency band- 
averaged SZ flux decrement in the frequency interval [86 GHz, 
94 GHz]. The ALMA telescope will observe at a frequency 90 
GHz with higher sensitivity and resolution than that of the GBT; 
comparison the sensitivity of the GBT with the ALMA telescope 
has already been studieqj. 

In Fig. [T] we show the dependence of the generalized spec- 
tral function G(x, T^) at a frequency 90 GHz (x= 1 .59), derived 
from Eq. (|5]l on temperature. We find that the absolute value of 
the generalized spectral function at a frequency 90 GHz signifi- 
cantly (and monotonically) decreases with electron temperature. 
The Kompaneets approximation is valid when the change of a 
photon frequency due to Compton scattering is much smaller 
than an initial photon frequency Av/v «: 1. Since the average 
photon frequency change per inverse Compton scattering equals 
Av/v = 4kbTJ(m^c^) (for non-relativistic electrons, see, e. g., 
Rybicki & Lightman 1979), the Kompaneets approximation is 
invalid when k^Te; ^ (meC-)/4, i.e. when kbTg ^ 128keV. Thus, 
strong deviations from the intensity value (and from the spectral 
function value of ^(1.59) ~ -3.28) derived in the Kompaneets 
approximation arise at high temperatures. The given estimate of 
the regime of breakdown of the Kompaneets approximation is an 
approximation and significant deviations (although not of order 
unity) take place at much lower temperatures, even at 30 keV 
(e.g. Fabbri 1981). 

An interesting feature of the SZ effect is that at frequency 
217 GHz (x=3.83) where the SZ effect in the framework of the 
Kompaneets approximation is zero, the SZ effect from an AGN 
cocoon is dominated by the high temperature electron compo- 
nent (T^e » 10*^ K). Therefore, the SZ effect at frequency 217 
GHz is an interesting tool for analyzing the hot electron compo- 
nent in AGN cocoons. Note that Colafrancesco (2005) has con- 
sidered a measurement of the SZ effect at a frequency 217 GHz 
to analyze a non-thermal electron component in radio lobes. 
Since the value of G(3.83, Te) is small when the temperature is 
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Fig. 1. Dependence of the generalized spectral function at fre- 
quency 90 GHz on the plasma temperature. 




http://safe.nrao.edu/wiki/bin/view/GBT/GBTSensitivityComparison 
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Fig. 2. Dependences of the generalized spectral function at the 
frequency 217 GHz and the combined generalized spectral func- 
tion C(2.26, 6.51, T) on plasma temperature are shown by solid 
and dashed lines, respectively. 

either much lower than 10** K or higher than 10'" K, we con- 
clude that this function should have an absolute maximum at an 
intermediate temperature. 

Dependence of the generalized spectral function G{x, T^) at 
a frequency 217 GHz (x=3.83) on temperature derived from Eq. 
(|5]l is shown in Fig. |2] The extremum of a curve G(3.83, Tg) is 
at a temperature ~ 160 keV (~ 1.9 x 10^ K) which is in the 
temperature range found from numerical simulations of cocoons. 
Therefore analysis of the SZ effect at this frequency should be 
a promising way of demonstrating the existence of an electron 
component at high temperature. Note that combined generalized 
functions defined in Sect. 2.2 are also shown in Fig.|2] 

Measurements at the frequency of 217 GHz are sensitive to 
the spectral response of the detectors. In Appendix B we show 
how the broad detector spectral response impacts on the possi- 
bility of an analysis of high temperature plasmas. 

In Sect. 2.2, we investigate the possibility of constraining 
the properties of high energy particle populations by means of 
SZ intensity measurements at two frequencies. 

2.2. Combined generalized spectral functions 

Measurements of the SZ effect at a frequency 217 GHz is a 
unique way of revealing a population of mildly relativistic elec- 
trons in AGN cocoons if the intensity distortion is observed at a 
single frequency. We propose below a new method for excluding 
the contribution from the low energy, non-relativistic electrons to 
the SZ effect by means of observations at two frequencies. 
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By C(xi,X2,Ti,), we denote a combined generalized spec- 
tral function which corresponds to a relativistic contribution to 
the generalized function at a frequency xi taking into account a 
measurement at a frequency X2 and is defined as 

CixuX2, T,) = G(xi, r,) - G(X2, T,) x (7) 

At low temperatures <K 10^ K, the limiting case (1) holds 
(see Sect. 2.1) and, therefore the combined generalized spectral 
function C{x\,X2, Tg) equals 0. It means that this function is in- 
sensitive to low-temperature electron components. (Note that a 
SZ intensity cavity at a frequency of 90 GHz produced by mildly 
relativistic electrons in an AGN cocoon due to relativistic SZ ef- 
fect corrections (Pfrommer et al. 2005) can vanish if the assump- 
tion of pressure equilibrium doesn't hold. This is because the SZ 
signal at a frequency of 90 GHz produced by a gas with a low 
temperature can be reproduced by a mildly relativistic gas with 
a high temperature and higher pressure, see Fig. [1]). Equation 
(|7]) defines a family of spectral functions: each choice of values 
for frequencies xi, X2 produces a different spectral function. The 
common property of spectral functions of this family is that they 
are more sensitive to mildly relativistic electron populations than 
to non-relativistic electron populations. Therefore, this choice 
for the combined generalized spectral function is motivated by 
studying of high energy electron populations. Let us consider 
different choices of frequencies of xi and X2. 

There are three basic spectral features that characterize the 
thermal, non-relativistic SZ effect signal: a minimum of its in- 
tensity located at a dimensionless frequency 2.26 (v=128 GHz), 
a crossover frequency xq = 3.83 (v=217 GHz), and a maximum 
of its intensity whose frequency location at a dimensionless fre- 
quency 6.51 (v - 369 GHz). If the xi - 3.83 then the values 
of generalized G(xi) and C{xi,X2) are coincident. If X2 - 3.83, 
the combined generalized function is undetermined because the 
value of g{x - 3.83) equals 0. We noticed that values of the 
spectral function of G{x, T^) at high temperatures at frequencies 
of xi - 2.26 and X2 - 6.51 cannot be simultaneously fitted by 
the spectral function of c x g{x), where c is the arbitrary con- 
stant. Therefore, the choice of frequencies xi = 2.26 and X2 - 
6.51 corresponding to minimum and maximum values of the SZ 
intensity in the Kompaneets approximation is suitable to analyze 
high energy electron populations. The dependence of the com- 
bined generalized function with xi - 2.26 and X2 = 6.51 on the 
plasma temperature is shown in Fig. |2] The combined function 
of C(1.59,6.51, Te) where the lower frequency corresponds to 
90 GHz is shown in Fig.|2]by dotted line. Figure |2] shows that 
the curves which correspond to the generalized spectral function 
at frequency 217 GHz and the combined generalized function 
C(2.26, 6.51, T) are very similar and have peaks at temperatures 
higher than 100 keV. Therefore, using the combined generalized 
function provides us with an alternative and equivalent method 
for studying a population of electrons with energies higher than 
100 keV. Such peaks at temperatures higher than 100 keV as 
that in the Fig. |2]permit us to maximize the SZ effect from very 
hot gas which is expected inside AGN cocoons from numerical 
simulations. 

To calculate the intensity map of the SZ effect using com- 
bined generalized spectral functions, we use Eq. (|4]i where the 
generalized spectral function G(x, Te) should be changed to the 
combined generalized spectral function C{x\,X2, Tg). The kine- 
matic Sunyaev-Zel'dovich effect which is a possible source of 
bias in the observations of the SZ effect by energetic electrons 
will be considered in Sect. 4. 1 . 



In the next section, an analysis of the SZ effect by means 
of the generalized spectral functions is considered for a specific 
astrophysical important case - the Sedov self-similar solution. 

3. A specific test case: the Sedov solution. 

Self-similar solutions for a strong point explosion in an ambient 
medium are used for modelling adiabatic supernova remnants, 
solar flares and processes in active galactic nuclei (Ostriker & 
McKee 1988). Sedov (1959) gives the analytical self-similar so- 
lution for description of the motion of a shock front and the 
distribution of fluid parameters inside the shocked region for a 
strong point explosion in a uniform ambient medium. The gas 
flow pattern is determined by only two parameters: the ambi- 
ent gas density pi, and the amount of energy E released in the 
explosion. The distance of the shock from the origin is given by 

where y6 is a numerical constant. 

The density p as a function of the dimensionless radial co- 
ordinate r/R decreases rapidly into the sphere and almost all the 
gas is in a relatively thin layer behind the shock wave (Sedov 
1959). As r/R — > 0, the pressure p tends to a constant limit and 
the temperature accordingly becomes very high. Since the SZ 
effect depends on the thermal energy density of the electron pop- 
ulation along a line of sight, a significant SZ effect is expected 
from the central region although the density is very small there. 

Yamada et al. (1999) investigated the SZ effect produced by 
cocoons of radio galaxies and constructed a model for the evolu- 
tion of a cocoon after the jet is turned off. They examined the 
evolution of a cocoon after the jet turns off by analogy with 
the evolution of a supernova remnant and showed that the co- 
coon remains hot enough to be the source of the SZ effect only 
during the free expansion and Sedov (adiabatic) phases. The 
Sedov stage of the evolution of the cocoon surrounding an AGN 
was also considered by Platania et al. (2002) and Chatterjee & 
Kosowsky (2007). In these papers, for simplicity the density of 
the gas inside the cocoon was assumed to be uniform and only 
the time evolution of the Comptonization parameter was taken 
into account. We will consider a more realistic case where the 
gas density depends on the position inside a cocoon. 
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Fig. 3. The intensity map |A///o| of the SZ effect derived from 
the assumption of the Sedov solution and the Kompaneets ap- 
proximation 

The values of pressure and temperature at any radius inside a 
cocoon are determined by the ambient gas density pi , the amount 
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of energy E released in the explosion, and time t. However, the 
parameters directly observed by X-ray telescopes are the am- 
bient number density n\ = pi/mp, the shock radius R and the 
Mach number M of the shock. The Mach number of a shock is 
usually derived from the Rankine-Hugoniot jump conditions (for 
a review, see Markevitch & Vikhlinin 2007). Another way to de- 
rive the Mach number of a shock is based on measurements of 
the flux ratio of the FeXXV and FeXXVI iron lines (Prokhorov 
2010). Since the amount of energy released and time can be ex- 
pressed in terms of the shock radius and the Mach number of 
the shock, we choose rii, R, and M as three parameters for de- 
scribing the Sedov solution. Typical values of ni = 10"~ cm"-', 
R - 5Q kpc, and M = 2.0 correspond to that of observed from 
AGN cocoons in galaxy clusters (Forman et al. 2005; Nulsen et 
al. 2005a). The Mach number of the shock can be calculated in 
terms of the shock velocity relative to the undisturbed gas given 
by (Sedov 1959) 



5 [pit^ 



1/5 



(9) 



The self-similar Sedov solution assumes that the pressure P2 
behind the shock is much larger than the pressure Pi of the am- 
bient gas (see Landau & Lifshitz 1959). The evolution of a su- 
pernova remnant in the high-pressure ambient gas was examined 
by Tang & Wang (2005) by means of high-resolution hydrody- 
namic simulations. They found that the supernova remnant evo- 
lution deviates from the Sedov solution when the shock radius is 
higher than the characteristic radius R^ (see Eq. 5 in their paper). 
Using the typical values of the number density «i - 0.01 cm"-', 
the ambient gas temperature T - 2 keV, and the release energy 
£ = 3 X 10^" erg, we found that the characteristic radius for 
an AGN cocoon equals 80 kpc and higher than the shock radius 
R - 5Q kpc adopted above. Therefore, the self-similar Sedov 
solution is a reasonable approximation in our case. 

Assuming that a short-lived AGN source can be modeled as 
an instantaneous central explosion, Cavaliere & Lapi (2006) cal- 
culated the SZ effect in the framework of the Kompaneets ap- 
proximation for the radial pressure profile taken form the Sedov 
solution. The intensity map of the SZ effect derived under the 
assumptions of the Sedov solution and the Kompaneets approx- 
imation is plotted in Fig. [3] We note that in this case, the inten- 
sity map of the SZ effect is almost plane with a slight intensity 
increase in the central region. 

To calculate the SZ intensity maps at frequencies 90 and 217 
GHz, we use the generalized spectral functions derived in the 
previous section for different values of temperature (see Eq. |5] 
Figs. 1 and 2). Using the Wright formalism and values of the 
three parameters introduced above, we calculate the SZ inten- 
sity map taking into account the presence of high temperature 
gas inside the cocoon. The intensity maps of the SZ effect at 
frequencies 90 GHz and 217 GHz are plotted in Figs. |4] and |5] 
respectively. 

Figure |4] shows the SZ intensity decrement at a frequency 
90 GHz at the center of this map in contrast to the SZ intensity 
map derived in the framework of the Kompaneets approximation 
(see Fig. [3]). This central SZ decrement is provided by electrons 
with very high temperatures presented in the central region of the 
cocoon and the existence of such a decrement is consistent with 
the decrease the generalized spectral function at a frequency 90 
GHz with temperature (see Fig.[T]). 

The SZ intensity map at a frequency 217 GHz shows a ring- 
like structure that is clearly seen in Fig. [5] This ring-like structure 
is associated with a SZ signal from a gas layer with temperature 
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Fig. 4. The intensity map |A///()| of the SZ effect at a frequency 
90 GHz derived from the Sedov solution in the framework of the 
Wright formalism. 



close to r a; 160 keV that corresponds to the maximum of the 
absolute value of the general spectral function (see Fig.|2]i. Since 
the value G(3.83, 160keV) » -1.2 of the generalized function at 
a frequency 217 GHz and temperature T a: 160 keV is much 
higher than that is at much lower and higher temperatures (see 
Fig. [2] and Eq. |6]), the SZ signal is therefore a promising tool 
for observing and analyzing plasma regions with a temperature 
close to 2 X 10'' K. 



3.1. Astrophysical application of tine ring at 21 7 GHz 

In this section, we consider a way of constraining the amount of 
energy E released during the active phase of a jet by using the 
location of a ring on the SZ intensity map at a frequency 217 
GHz (or on the SZ intensity map derived from measurements at 
the 128 and 369 GHz frequencies, see Sect. 2.2). 




X (kpc) 



Fig. 5. The intensity map |A///o| of the SZ effect at a frequency 
217 GHz derived from the Sedov solution in the framework of 
the Wright formalism. 



The gas temperature at the back of the shock can be ex- 
pressed in terms of the shock velocity u as (Sedov 1959) 



2{y — l)mpM^ 

hiy + 1)2 



(10) 



where nip is the proton mass and y is the adiabatic index. The 
time-dependence of the gas temperature behind the shock is then 



T2 = 



- l)mpE^'' 
25{y + l)^hp]'^t^/^ ' 



(11) 
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The temperature profile inside the sphere (Landau & Lifshitz 
1959) is written as 



7^^3/(7-1) 



(12) 



where r is the distance from the center Although Eq. (fT2l i was 
derived under the assumption that r «: we verify that this 
equation can also be applied if r < R. Since the maximal con- 
tribution to the SZ effect at a frequency 217 GHz comes from 
the plasma with temperature close to k\,Ti - 160 keV, we fix this 
temperature and derive the dependence of the location of a ring 
on time 



'f2(2-y)/5£(2r+l)/15' 



r 

where 



■(r~i)/3 



Pi 



(2r+i)/i5 



A = 
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If the adiabatic index equals y 
given by 



n2/9 



2 £■13/3 



13/3 



1/15 



5/3 then the radius of a ring is 



(15) 



Note that the ring radius depends weakly on time r oc r^^'^ 
and represents a parameter which changes negligibly with time 
(in contrast to the shock radius, see Eq. |8]i. The radius of a ring 
can be used as one of the parameters for describing the Sedov 
solution. If the ring radius is used instead of the Mach number 
then the value of the released energy E is derived from Eqs. ([8]|, 
( fTSl ), depends strongly on the ring radius, and is given by 



E oc 



RV2 ■ 



(16) 



Since the ring radius is r{t) oc f^/'^ and the distance of the shock 
from the origin is R{t) oc f-/^, the ratio of P{t) to R{t) is constant 
in time and the value of the energy is essentially determined by 
the ring radius of the SZ signal at a frequency 217 GHz. 

The analytical solution of hydrodynamic problems such as 
the Sedov solution is useful, but many problems have no such 
solution and numerical simulations are needed. Furthermore, the 
moderate values observed for Mach numbers of shocks require 
inclusion of the effects of gravity and finite initial gas pressure in 
the pre-shock regions. We consider the SZ effect by relativistic 
electrons residing in AGN cocoons during the active jet phase 
by using the results of numerical simulations in Sect. 4. 

3.2. Synthetic observations of tlie SZ effect on AGN cocoons 
in gaiaxy dusters 

High-resolution observations of the SZ effect are a promising 
tool for observing AGN cocoons in galaxy clusters. Powerful 
AGN outbursts were revealed in the Hydra A cluster (Nulsen et 
al. 2005b), in the Hercules A cluster (Nulsen et al. 2005a), and 
in the MS0735.6-H7421 cluster (McNamara et al. 2005). These 
clusters are located at the redshifts of z=0.05, 0.15 and 0.22, 
respectively. A flat ACDM cosmology, with Hq = 70 km s"' 
Mpc"^ and Qm - 0.27, gives a scale of 2.67 kpc arcsec"' for 
the redshift of z=0.15. Therefore, with the spatial resolution of 



ALMA (2.8" at 1 10 GHz), we will be able to resolve the struc- 
tures down to 10 kpc at z=0.15. In this section we simulate ob- 
servations of the ring-like structure on the SZ intensity map stud- 
ied in the previous section. 

An ALMA observation of the central region of the Perseus 
cluster is simulated by Pfrommer et al. (2005) to reveal a pres- 
ence of cavities on the SZ intensity map. We convolve the SZ 
intensity map with a Gaussian to obtain the resolution of the 
ALMA compact configuration as it was done by Pfrommer et al. 
(2005). Since the angular extent of a 100 kpc box (considered 
in Sect. 3) at the redshift z=0.15 is close to the ALMA field of 
view (^ 30"), for the sake of simplicity we put a synthetic AGN 
cocoon in a galaxy cluster at the redshift z=0.15. The resulting 
synthetic SZ intensity map of an AGN cocoon at a frequency of 
217 GHz in a galaxy cluster at z=0.15 is shown in Fig.|6l Figure 
|6] shows that the ring-like structure is still present on the SZ in- 
tensity map after a convolution with Gaussian of FWHM^ 3.0" 
and, therefore, the ring radius measurement will be an interesting 
method provided by the next generation telescopes (e.g. ALMA) 
to determine the energy released during an outburst. The spatial 
resolution of 3" is sufficient to resolve the ring-like structure on 
the SZ intensity map at frequency of 217 GHz and, therefore, 
ALMA with the spatial resolution of 1.3" at a frequency of 217 
GHz in the most compact configuration should reveal the ring- 
like structures on SZ intensity maps. If we choose the same AGN 
cocoon in a cluster at z=0.05, the spatial resolution of 10" is suf- 
ficient to reveal the ring-like structure on the SZ intensity map. 




X (kpc) 



Fig. 6. Synthetic observation of an AGN outburst in a galaxy 
cluster The intensity map |A///o| of the SZ effect at a fre- 
quency 217 GHz derived from the Sedov solution in the frame- 
work of the Wright formalism and smoothed to the resolution of 
FWHM^ 3.0". 



We simulated SZ maps for ALMA with the spa- 
tial resolutions of 3" and 1.3" at frequencies of 90 
GHz and 217 GHz, respectively. Our simulated SZ in- 
tensity maps for ALMA are similar to those in Figs, 
m and ID and available for download at the webpage 
|http : //www ■ pact . inaf . it/cosmoct7web_group/research47cocoon_sz . html| 

Pfrommer et al. (2005) examined whether the plasma bub- 
bles are detectable by ALMA and predicted 5 cr detection of the 
bubble of the Perseus cluster in an exposure time of 5 hours. 
Note that a signal at a frequency of 217 GHz derived from 
the Sedov solution is four times smaller than that analyzed by 
Pfrommer et al. (2005) at a frequency of 144 GHz (see Figs. 
|5]of this paper and 1 of Pfrommer et al. 2005). The signal-to- 
noise ratio for a detection of an AGN cocoon is proportional the 
square root of the exposure time (e.g. Pfrommer et al. 2005; for 
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a review, see a document 0) and, therefore, to detect the signal 
at a frequency of 217 GHz the ALMA exposure time of 80 (i.e. 
5 X 4^) hours is required. The ALMA Design Reference Science 
Plan (Version 2.2) contains two projectfl^in SZ studies (1.4.1 
and 1 .4.2) with integration time for each program of 400 hours 
and new projects concern the SZ effect on AGN cocoons should 
be promising. The Chandra X-ray satellite with a high spatial 
resolution of 0.5" has studied a non-relativistic electron compo- 
nent in the ICM, while ALMA high-resolution observations of 
AGN cocoons in galaxy clusters should permit us to reveal a new 
(mildly relativistic) electron component. 



4. The SZ effect from a simulated cocoon 

Recent hydrodynamic simulations do not take into account rel- 
ativistic corrections to the SZ effect to derive SZ intensity maps 
(e.g. Sijacki et al. 2008). Such relativistic corrections are neces- 
sary for calculating the SZ effect on AGN cocoons produced by 
electrons with high temperatures. So far, studies of the SZ effect 
in the relativistic formalism were done only for an analytic toy 
cocoon model (see Pfrommer et al. 2005). In this section, we de- 
rive the SZ intensity map by using hydrodynamical simulations 
and the Wright relativistic formalism. 

In previous papers (Antonuccio-Delogu & Silk 2008; Tortora 
et al. 2009), we used an Adaptive Mesh Refinement (AMR) code 
to follow the evolution of the cocoon produced by the jet prop- 
agating in the interstellar (ISM)/intergalactic medium (IGM). 
To perform the simulation, we used FLASH v.2.5 (Fryxell et 
al. 2000), a parallel, AMR code, which implements a sec- 
ond order, shock-capturing, Piecewise-Parabolic-Method (PPM) 
solver The jet is modelled as a one-component fluid, character- 
ized by a density which is a small fraction of the initial density 
of the ISM. In this simulation, the power injected by the jet is 
10"*^ erg/s. We model the environment, where the jet propagates, 
as a two-phase ISM, comprising a hot, diffuse, low-density com- 
ponent having a temperature and a central density at 10^ K and 
1 cm^, respectively, and a cold, clumped system of high-density 
clouds in pressure equilibrium with the diffuse component. Such 
values of temperature and pressure are typical of the ISM in the 
central parts of an elliptical at high redshift. We assume that the 
diffuse gas is embedded within a dark matter (DM) halo, the lat- 
ter being described by a Navarro-Frenk- White (NFW) density 
profile. 

We choose a simulation box having a size Lbox = 40/!_i kpc, 
where h is defined such that the Hubble constant Ho is 100 h km 
s ' Mpc The spatial resolution attained is a function of the 
maximum refinement level and of the structure of the code. For 
a block-structured AMR code like FLASH, where each block is 
composed by n ^ x cells, the maximum resolution along each 
direction is given by Lbox/(«A2'), where I is the maximum refine- 
ment level. In this simulation, n ^ - = 8 and I - 6, thus the 
minimum resolved scale is 78. 125 pc. Note that we are perform- 
ing a 2D simulation, but we do not impose any special symmetry. 

Numerical simulations by Antonuccio-Delogu & Silk (2008) 
and Tortora et al. (2009) show that the temperatures can reach 
very high values within the cocoons (T x 10^ - 10" K), if the 
ISM is dominated by a population of cold, star forming clouds 
embedded into and in approximate pressure equilibrium with a 
hot, diffuse phase. The pressure within the cocoon can reach high 
values because the temperatures are on average very high. The 



SZ effect should be significant since it determined by the pres- 
sure of the electron population integrated along a line of sight. 
Our simulation temperature and pressure maps are shown in Fig. 
HJandl 



20 
X (kpc) 



Fig. 7. The simulation map of the gas temperature (in K) in the 
AGN cocoon 



Figure Q reveals a hot shell with a temperature of k^Tg > 2 
MeV around the internal region of the simulation AGN cocoon. 
The temperatures of the internal region are smaller than that of 
the hot shell and, therefore, absolute values of the generalized 
spectral functions at frequencies of 90 and 217 GHz should be 
higher for the internal region in the corresponding temperature 
range (see Figs. [1] [2] and|7]i. 




' http://www.cv.nrao.edu/naasc/ALMAsensitivity.ps 
http://www.eso.org/sci/facilities/alma/science/drsp/ 



10 20 30 40 

X (kpc) 

Fig. 8. The simulation map of the gas pressure (in erg/cm^*) in 
the AGN cocoon 

We find that the pressure values are higher inside the inter- 
nal region in comparison with the thermal hot shell pressure. 
Therefore, the SZ effect in the cocoon internal region should 
be more significant than that produced by electrons located in- 
side the hot shell. The pressure inside the simulation cocoon is 
two orders of magnitude higher than that found in the diffuse 
ISM. Such a pressure-jump corresponds to the Mach number of a 
shock equaled M ^ yfAPTfWx ~ 10 (Landau & Lifshitz 1959). 
The plasma inside the AGN cocoon is largely over-pressured 
relative to the ambient ISM. The radio lobes of the nearby ra- 
dio galaxy Centaurus A is expanding into the ISM at the Mach 
number of 8.5 (Kraft et al. 2003) is representative of such over- 
pressured plasmas. 

To produce the 3D pressure and temperature maps which are 
necessary to derive the intensity map of the SZ effect we rotate 
the 2D pressure and temperature maps along the jet axis. We cal- 
culated the SZ effect using the values of the generalized spectral 
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functions found in Sect. 2. Since the gas temperatures in a co- 
coon are high, such a cocoon should be a source of the SZ effect 
at a frequency 217 GHz. The intensity map of the SZ effect at a 
frequency 217 GHz derived from the simulation maps of the gas 
pressure and temperature is plotted in Fig.|9] 



20 

(kpc) 



Fig. 9. The intensity map |A///o| of the SZ effect at a frequency 
217 GHz derived from the numerical simulation in the frame- 
work of the Wright formalism 



A comparison of Figs.|2]and[8]with Fig. |9] strongly suggests 
that the dominant contribution to the SZ effect at a frequency 
217 GHz originates from the cocoon internal region. Therefore, 
measurements of the SZ effect at a frequency 217 GHz are a 
powerful tool for potentially revealing the dynamically domi- 
nant component inside AGN cocoons (see Fig.|8]l. The proposed 
method based on the generalized spectral function is a way of 
maximizing the contribution from the gas with a temperature in 
the range 10'' K < Te < 10'" K (see Fig. |2|. The study of the 
SZ effect in a plasma at such temperatures is important since 
numerical simulations predict that such a plasma represents the 
still invisible dynamically dominant component inside AGN co- 
coons. 

Our simulations are characterized by ISM density and tem- 
perature at 1 cm"-' and 10^ K (ISM pressure at ~ 10"^ erg/cm-'), 
respectively, typical of the ISM in the central parts of an ellip- 
tical at high redshift (z ^ 1). We found that the plasma inside 
the AGN cocoon is largely over-pressured relative to the ambient 
ISM. The derived intensity A/ of the SZ effect in this case is three 
orders of magnitude higher than that of cocoons in galaxy clus- 
ters. Though the linear scales of cocoons are close in both cases, 
the angular diameter of the cocoon within the distant elliptical 
galaxy equaled to 0z=i = Lbox/^^A ~ 2.5" is much smaller than 
that of the cocoon with the same linear size inside the Perseus 
cluster 0z=o.oi79 = 1-8', the angular distance is denoted by Da- 
Then the SZ fluxes from such cocoons which are proportional to 
AI/D^ should be of the same order (see Figs.|5]and|9l). 

The main difference in the intensity of the signal is then a 
result of the higher internal pressure in the simulated cocoon. 
The pressure of high temperature gas in the cavity center, in the 
framework of the analytical model in Sect. 3, can not be much 
higher than the ambient pressure of a galaxy cluster ~ 3 x 10"" 
erg/cm"' because of the moderate Mach number of the shock of 
M -2 and a decrease of the pressure from the shock front to the 
cavity center (e.g., Sedov 1959). 



4.1. An analysis of the hot gas by means of the combined 
generaiized function 

As was shown in Sect. 2.2 using the combined generalized spec- 
tral functions provides us with an alternative to a measurement 
of the SZ effect at a frequency 217 GHz for studying a pop- 
ulation of electrons with energies higher than 100 keV. In this 
section we derive the intensity map of the SZ effect by means 
of the combined function C(2.26, 6.51, Tg) and consider how the 
contribution from the kinematic Sunyaev-Zel'dovich (SZ) effect 
can be eliminated (for review of the kinematical SZ effect, see 
Birkinshaw 1999). 

The generalized spectral function at a frequency 217 GHz 
belongs to the family of spectral functions defined by Eq. 
Q. Both G(x, Tf.) and C(2.26,6.51, T^) functions have absolute 
maximal values at temperatures above 100 keV and, therefore, 
these functions are interesting for analyzing the presence of high 
energy electrons which are expected to be inside AGN cocoons. 
The intensity map AI/Iq of the SZ effect derived from our nu- 
merical simulation by using the combined generalized spectral 
function C(2.26, 6.51, T) is plotted in Fig.fTUl 




X (kpc) 



Fig. 10. The intensity map |A///o| of the SZ effect derived from 
our numerical simulation by using the combined generalized 
spectral function C(2.26, 6.5 1 , T) 

The similarity in shape between the SZ intensity maps in 
Figs.l9land[T0lshows that a population of high energy electrons 
can be constrained by means of different and equivalent methods 
based on measurements of the SZ effect. 

A source of bias in the observations of the SZ effect by en- 
ergetic electrons could be provided by a possibly relevant kine- 
matic SZ effect. The kinematic SZ effect arises from the bulk 
motion of the medium relative to the CMB rest frame. The 
change of the CMB intensity due to the kinematic SZ effect is 
given by 



Alidnix) = -PTh(x), 



(17) 



where y6 is the line-of-sight velocity of bulk motion in units of 
the speed of light. The spectral function h{x) is given by 



h(x) = 



x^ exp(x) 
[exp(x)-l]2- 



(18) 



The absolute maximum of the kinematic SZ effect lies at 
a frequency 217 GHz which coincides with the crossover fre- 
quency of the non-relativistic thermal SZ effect. Pfrommer et 
al. (2005) concluded that it is impossible to remove the degen- 
eracy of the cocoon optical depth and the kinematic SZ effect 
using only a single frequency SZ observation. Colafrancesco 
et al. (2009) found that the value of the slope S - (_AI{x^) - 
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A/(xb))/(xa - jcb) of the SZ effect does not depend on the kine- 
matic SZ effect spectrum in the frequency range around the 
crossover frequency of the thermal SZ effect, i.e. in the fre- 
quency range x^,x\, 6 (3.5,4.5). We now show that there is 
another way to eUminate the kinematic SZ effect based on the 
combined generaUzed spectral functions. 
Let us define the intensity function U as 



g(-^min) 



(19) 



A/(xo) X 



n) 



h(xo) Kxq) g(jCmax) 



where AI{x) is the total intensity of the thermal and kinematic SZ 
effects, jc„in = 2.26 (v = 128 GHz), xq = 3.83 (v=217 GHz), and 
Xm-dx = 6.51 (v = 369 GHz) correspond to the frequencies of the 
minimum, the crossover and the maximum of the non-relativistic 
thermal SZ effect, respectively. 

Since A/(Xmin) - AI{x^^^) x gixminJ/gixmax) docs not depend 
on the non-relativistic thermal SZ effect (see Sect. 2.2) and ^(.yo) 
equals zero, there is no contribution of the non-relativistic ther- 
mal SZ effect to the intensity function U. Note that AI(x„iin) - 
AI{xq) X h{x^i„)/li(xo) and AI{x^^^) - AI{xq) x li(x^dx)/ Kxq) 
have no dependence on the kinematic SZ effect, and therefore 
the study of the intensity function [/ is a way of constraining 
a population of high energy electrons without biases provided 
by a possible relevant kinematic SZ effect. The value of the last 
bracket on the right-hand side of Eq. ( fT9] l is very close to 1 (it 
equals 1.029) and, therefore, the spectral function GU{T) which 
corresponds to the intensity function U (see Eq. |4]i is approxi- 
mately equal to 



GU{T^) ~ C(2.26,6.51,re)-G(3.83,re) 



(20) 



where the spectral functions C(2.26, 6.51, Te) and G(3.83,re) 
correspond to the SZ intensities of A/(;iCniin) - A/(x,nax) x 
gixmin) / gixmax) and AI(xo), respectively. The spectral function 
GU{T) is shown in Fig. [TT] 



In the next section, we consider the scattered CMB photon 
spectrum at higher frequencies. This study is important since 
photon energy gains due to the inverse Compton effect on high 
energy electrons are substantial. 

5. Flatness of the high-frequency tail of the 
scattered photon spectrum 

Electrons with high energies produce a substantial energy gain 
of photons, which scattered to the high-frequency tail. For an 
optical depth r <*: 1 and a sufficiently large Comptonization pa- 
rameter y, a sequence of declining peaks in the high-frequency 
tail (jc » 10) of the scattered photon energy flux spectrum is 
present (e.g., see Loeb et al. 1991). These peaks correspond suc- 
cessively to singly, doubly (etc) scattered photons (see Figs. 4a- 
4d in Loeb et al. 1991). Next, we study the high-frequency tail 
of the scattered CMB spectrum by high energy electrons at fre- 
quencies below the first peak frequency. 

At high frequencies (i.e. x » 10) the term x^/(exp(x) - 1) 
of Eq. (|5]) decreases strongly and, therefore, the Eq. (|5]) can be 
written as 



G(X, Te) 



f 



Pi(s,T,) x^exp(-3^) 
©(Te) exp(xexp(-i)) - 1 



-ds. 



(21) 



The treatment is considerably simplified if exp(xexp(-i)) » 1 
and in the case the generalized spectral function G{x, Te) is given 
by 



G{x, T,) 



f 



x^ exp (-3s - xe\p(-s)) ds. 



(22) 



The sub-exponential function f{s) - -3s - xexp(-s) has 
a maximum at the frequency shift s^ax = ln(A/3). Since 
exp(xexp(-Sn,ax)) » 1 the approximate expression for the gen- 
eralized spectral function is valid and we calculate the integral in 
Eq. (l22l) by the Laplace's method. The approximate value of the 
sub-exponential function in a neighborhood of the points = s^^x 
is then 




Fig. 11. Dependence of the function GU{Te) at on the plasma 
temperature. 

The intensity map Al/ Iq of the SZ effect derived from our nu- 
merical simulation by using the function GC/(r) corresponds to a 
residual between the intensity maps of Figs.[T0]and|9l The func- 
tion GU{T) has the minimum at temperatures higher than 300 
keV and is not affected by the kinematic SZ effect, and therefore 
differencing these SZ intensity maps (see Fig. l9l andfTOb is an 
important tool for analyzing a high energy electron population. 



f{s) = -3Sn,ax - Xe - -e X (S - Sn,ax) 

and the spectral function approximately equals 

i(lnf,re) 
©(T'e) ■ 



(23) 



(24) 



As was shown by Kino et al. (2007) and Antonuccio-Delogu 
& Silk (2008), high gas temperatures (^bT'e ~ 1 MeV) are ex- 
pected in AGN cocoons. For example, we choose the tempera- 
ture equaled to 500 keV to study the high-frequency tail of the 
generalized spectral function. For high temperatures the function 
Pi{s, Tg) is wide and it is centered at high values of the frequency 
shift s. The distribution of frequency shifts for single scattering 
Pi{s, Ts) at the temperature of kbTf. = 500 keV is shown in Fig. 

Figure [12] shows that the values of the distribution of fre- 
quency shift lie in the narrow range (0.25, 0.28) if the frequency 
shift s is in the range (1.5, 3.0). The frequency shift s^ax lies in 
this range when the dimensionless frequency is of 15 < x < 90 
and, therefore, the approximate value of the generalized spectral 
function in this frequency range is 



G(x,500keV) ~ 0.25 



— (1 

T\7 



1 



©(500keV) 



(25) 
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Fig. 12. Dependence of the distribution of frequency shifts 
Pi{s,k\,Tf.) at the temperature of 500 keV. 

The quantitative dependence of the spectral function 
G(x, khTfi = 500 keV) on the dimensionless frequency x is il- 
lustrated in Fig. [13j 
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Fig. 13. Dependence of the spectral function G{x,k-oTs - 500 
keV) on the dimensionless frequency x is shown by the solid 
line. The approximate value of G (x, fcbTe - 500 keV) is shown 
by the dashed line. 

We conclude that the generalized spectral function should be 
flat in the broad frequency range when the temperature values 
are sufficiently high. Since the non-relativistic spectral function 
g{x) is a rapidly decreasing function at high frequencies {x > 10) 
in contrast with the spectral function G(x, T^) at high electron 
temperatures, a measurement of the SZ efifect at high frequen- 
cies provides an interesting test of the presence of high energy 
electrons. 

6. Discussion and conclusions 

The direct detection of the SZ eff'ect can in principle provide 
a unique diagnostic tool to study the physical conditions of 
the ICM at high redshift. To reach this end, an accurate treat- 
ment of the spectral properties of the SZ signal is very impor- 
tant, and this was our main aim in this work. Previous models 
were based on simplifying assumptions, like the assumption of 
pressure equilibrium and density homogeneity (Colafrancesco, 
2005; Pfrommer et al., 2005). Here we have instead considered 
the SZ signal arising from a inhomogeneous exactly solvable 
configuration, i.e. a spherically symmetric Sedov-expanding re- 
gion, and even more realistic models obtained from 2D fluid- 
dynamical simulations of jet/cocoon system propagating into the 



ISM/IGM. 

In the standard model of a typical strong double radio sources, 
a couple of twin jets advances into the IGM and engenders a 
hot, low-density "cocoon" which is filled only of matter originat- 
ing from the jet and of part of the shocked IGM (Scheuer 1974; 
Blandford & Rees 1974). Numerical simulations by Antonuccio- 
Delogu & Silk (2008), Tortora et al. (2009) and theoretical works 
by Kino et al. (2007) suggest that the plasma temperatures within 
the cocoons are very high {T ^ 10^ - 10" K) and, therefore, 
these plasmas are invisible in the soft X-ray band and generate 
cavities on the X-ray maps (e.g. McNamara et al. 2005). 

The energy density of the cocoon is higher than that of the 
IGM: it is a local hot region that induces a significant thermal 
SZ effect (Yamada et al. 1999). Since the temperature within the 
cocoons is expected to be very high, we have adopted a rela- 
tivistically corrected formalism to calculate the SZ eff'ect. Using 
the relativistically corrected Wright formalism, Colafrancesco 
(2005) and Pfrommer et al. (2005) have shown that the hot 
plasma within the X-ray cavities should be observable by detect- 
ing the SZ signal at a frequency v =217 GHz, since the contribu- 
tion of the non-relativistic IGM to the SZ effect at this frequency 
is negligible. 

However, as stated above, the analytical models for the co- 
coons considered by Colafrancesco (2005) and Pfrommer et al. 
(2005) are based on the assumptions of pressure equilibrium 
and cocoon homogeneity. On the other hand, recent numerical 
simulations of the SZ effect do not include the relativistically 
corrected SZ treatment (e.g. Sijacki et al. 2008). To produce 
more realistic SZ maps we have considered the analytical Sedov 
model and numerical 2-D simulations and calculate the SZ effect 
by using the Wright formalism. 

We have studied the CMB distortion due to the SZ effect as 
a function of the gas temperature. We have shown that the CMB 
spectral distortion can be expressed in terms of the generalized 
spectral functions G{x, 7^) which are functions of the gas tem- 
perature. This spectral function at the frequency 217 GHz has a 
peak at a temperature kbTf. - 160 keV and, therefore, analysis 
of the SZ effect at this frequency should be a powerful way of 
checking for the existence of electron populations at high ener- 
gies. 

We have also proposed a new method for excluding the contri- 
bution from the low energy, non-relativistic electrons to the SZ 
eff'ect by means of observations at two frequencies. The derived 
similarity of the generalized spectral function G(3.83, Te) and of 
the combined generalized function C(3.83, 6.51, Tg) provides us 
with an alternative and equivalent method for studying a popu- 
lation of electrons with energies higher than 100 keV. We have 
shown how one may correct for a possible contribution from the 
kinematic SZ eff'ect using the difference between the general- 
ized spectral function G(3.83, Tg) and the combined generalized 
function C(3.83, 6.51, Tg). 

In Sect. 3 we have studied in detail the SZ effect for a spe- 
cific astrophysical important case - the Sedov self-similar solu- 
tion. The Sedov solution was used by McNamara et al. (2005) 
to described the observed X-ray cavities. We have then consid- 
ered the AGN cocoon in a galaxy cluster, and calculated the SZ 
effect at a frequency 217 GHz produced by the hot gas inside an 
X-ray cavity. We have used the radius of the ring structure on the 
SZ intensity map for deriving the value of the released energy. 
In Sect. 3.2 we have also showed that the ALMA telescope is 
needed to reveal mildly relativistic electron populations in AGN 
cocoons. 

We have then proceeded by considering the SZ effect from a 
simulated cocoon in Sect. 4. The simulated cocoon is largely 
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over-pressured relative to the ambient ISM. We have derived 
the SZ intensity map by using hydrodynamical simulations and 
the Wright relativistic formalism. We showed that the SZ fluxes 
from cocoons in the central part of a distant elliptical and a 
nearby galaxy cluster are of the same order. 

We have showed that the CMB spectral distortion is flat in 
a broad range of high frequencies if gas temperatures are high, 
and we have estimated the approximate value of the generalized 
spectral function at temperature 500 keV in the frequency range 
15 < X < 90. 

Our study demonstrates that the SZ effect, originating from 
populations of high energy electrons residing in AGN cocoons, 
calculated using the SZ intensity maps in the framework of rel- 
ativistically corrected formalism, can give a detectable signal. 
We have proposed methods for maximizing the contribution of 
high energy electrons, which provides promising tests for find- 
ing the hitherto invisible hot component of AGN cocoons. The 
detections and measurements of these populations are within the 
reach of ALMA, and could provide a fundamental tool to char- 
acterize the evolution of AGN activity at high redshifts. 



Appendix A: Fluid dynamical scaling of cocoon 

In the Appendix A we show how one can produce intensity maps 
of the SZ effect by means of a dynamical scaling of cocoon. We 
will use a slightly modified version of the scaling fluid dynami- 
cal equations by Tang & Wang (2009), which applies to the case 
of an expanding cocoon continuously powered by a jet. The set 
of equations we solve is very similar to Eqs. (l)-(4) of Tang & 
Wang (2009), with the only difference that the energy equation 
now contains a source term describing the injection of energy 
density: 



dpe 



H- V ■ [(pe + P) v] = Jd{r -n)-L 



(A.l) 



where: J is the power injected by the jet into the cocoon, S{r-rh) 
denotes a 3-dimensional Dirac's delta function with ry, denoting 
the (bounded) jet's injection region, and L = iigiiHMT) is the 
cooling term. Tang and Wang (2009) assume a scaling of the 
form: q2 - qiQ''', where: q - (L, M, t, T, v, e,p, p, E, J) is either 
a fundamental or derived physical quantity: length L, mass M, 
time f, temperature T, velocity v, specific energy e, density p, 
pressure p, total energy E, and jet's power J) and Q is a positive 
constant. Under the hypothesis that the cooling function A ^ 
r'/-, which holds for T > 5 x 10'' K in a fully ionized plasma, 
the scaling quantities are found to be connected by the following 
relations: 



Eq. IA.8I was derived from the scaling of the cooling term, and 
the final equation (eq. lA.9b derives assuming the following scal- 
ing of jet's power: J oc JoQ'' . 

Note that Eq. IA.8I is different from Eq. (15) of Tang & Wang 
(2009). The latter reads: im - -3it + 5/l + 0.5/7-. Our derivation 
is obtained starting from the scaling of the right-hand side terms: 
d (pe) Idt ~ ip + ie - it - im - 3/l + - i,) - i, - Im - it - 3/,, 
where the symbol "~" denotes the power law scaling index. The 
number density scales as « ~ p and, therefore, the cooling term 
scales as L ~ 2/m - 6//, + Q.Sir- Thus, by equating the two terms, 
we have: im - ii - 3// = 2;^ - 611 + 0.5;V, which gives Eq. IA.8l 
As described by Tang & Wang (2009), this linear system is com- 
pletely specified after the choice of 2 of the coefficients. We 
choose to fix the scaling coefficients for length and time, ii and 
and find 



V = 'L - 


h 


e = 'T = 


2/l - 2 


p = iL 


2i, 


- = 3/l 


-4/t 


E = 6/l 


-4/, 


M = 4/l 


-2/, 


■j = 6/l - 


- 5/,. 



(A. 10) 



Thus, a given solution (also numerical) can be rescaled by 
means of Eqs. dA.lOb to a different jet's power using scalings 
coefficients. 

To produce an intensity map of the SZ effect by a cocoon 
surrounding a less powerful jet using the dynamical scaling of 
cocoon we choose scaling parameters equaled: Q = 10, //. = 
-0.2, and /, = -0.05. The scaled power injected by the jet is lO'*^ 
erg/s and the scaled size of the simulation box is 25/i_i kpc. The 
intensity map of the SZ effect at a frequency 217 GHz derived 
from the scaling of the simulation maps of the gas pressure and 
temperature is plotted in Fig. lA.ll 



h — Ih — h 



(A.2) 



/e - 2/l - 2/t 

ip - iu - 3/l 
ip - iu - ih - 2/t 
/e = /m + 2/l — 2/t 
/t = 2/l - 2/t = /e 
/m = -3/t + 5/l - 0.5/t 
/j — /m + 2/l — 3/t 



(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 



(kpc) 



Fig. A.l. The intensity map |A///o| of the SZ effect at a frequency 
217 GHz derived from the dynamical scaling. 



Thus a new intensity map of the SZ effect can be found for 
cocoons evolving in different ISM by means of one particular 
simulation and of the dynamical scaling. Therefore, the dynam- 
ical scaling is a promising method to produce the intensity maps 
of the SZ effect by hot gas residing in AGN cocoons in different 
ISM. 
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Appendix B: Measurements at 217 GHz 

In Appendix B we show how the broad detector spectral re- 
sponse impacts on the possibility of an analysis of high tem- 
perature plasmas. 

We convolve the spectral function g{x) over a Gaussian de- 
tector response centered at 217 GHz as it was done by Shimon 
& Rephaeli (2003). The signal at a frequency of 217 GHz from 
a low-temperature plasma which arises due to the convolution 
over a spectral response is shown in Fig. B. 1 as a function of the 
FWHM of the Gaussian detector response. 



0.15 




0.00 r I 

10 20 30 '10 

FWHM (GHz) 



Fig.B.l. The contribution of the SZ signal at a frequency of 217 
GHz due to the convolution over a Gaussian spectral response as 
a function of the FWHM 
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The contribution of the SZ signal at a frequency of 217 GHz 
due to the convolution over a spectral response with the FWHM 
< 30 GHz is far smaller than the SZ signal due to the presence 
of high temperature plasmas (see Fig. |2]i and, therefore, the SZ 
effect at frequency of 217 GHz is a promising tool for analyzing 
the hot electron component in an AGN cocoon if the FWHM of 
a spectral response is smaller than 30 GHz. 
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